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Abstract 

In this paper, we analyze the symbol error rate (SER) of space-time network coding (STNC) in 
a distributed cooperative network over independent but not necessarily identically distributed (i.n.i.d.) 
Nakagami-m fading channels. In this network, multiple sources communicate with a single destination 
with the assistance of multiple decode-and-forward (DF) relays. We first derive new exact closed-form 
expressions for the SER with M-ary phase shift-keying modulation (M-PSK) and M-ary quadrature 
amplitude modulation (A/-QAM). We then derive new compact expressions for the asymptotic SER to offer 
valuable insights into the network behavior in the high signal-to-noise ratio (SNR) regime. Importantly, we 
demonstrate that STNC guarantees full diversity order, which is determined by the Nakagami-m fading 
parameters of all the channels but independent of the number of sources. Based on the new expressions, 
we examine the impact of the number of relays, relay location, Nakagami-m fading parameters, power 
allocation, and nonorthogonal codes on the SER. 
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I. Introduction 

Cooperative communication has been recognized as a promising low-cost solution to combat 
fading and to extend coverage in wireless networks [1,2]. The key idea behind this solution is to 
employ relays to receive and transmit the source information to the destination, which generates a 
virtual multiple-input and multiple-output (MIMO) system to provide spatial cooperative diversity 
[3,4]. Apart from diversity, throughput enhancement is another critical challenge for wireless 
networks. Against this background, network coding (NC) is proposed as a potentially powerful 
tool to enable efficient information transmission, where data flows coming from multiple sources 
are combined to increase throughput [5-7]. 

Recently, the joint exploitation of cooperative diversity and NC has become a primary design 
concern in distributed networks with multiple users and multiple relays [8-11]. Motivated by 
this, [8] investigated various sink network decoding approaches for the network with intermediate 
node encoding. In [9], linear network coding (LNC) was applied in distributed uplink networks to 
facilitate the transmission of independent information from multiple users to a common base station. 
In [10], low-density parity-check (LDPC) code and NC was jointly designed for a multi-source 
single-relay FDMA system over uniform phase-fading Gaussian channels. In [11], cooperative 
network coding strategies were proposed for a relay-aided two-source two-destination wireless 
network with a backhaul connection between the sources. 

In order to increase the capacity and transmission reliability of wireless cooperative networks, 
multiple antennas are deployed to gain the merits of MIMO processing techniques [12-17]. In this 
strategy, distributed space-time coding (DSTC) was proposed to further boost network performance, 
where the antennas at the distributed relays are utilized as transmit antennas to generate a space-time 
code for the receiver [18, 19]. A differential DSTC was proposed in [20] to eliminate the requirement 
of channel information at the relays and the receiver. The combined benefits of maximum-ratio 
combining (MRC) and DSTC were investigated in [21]. In [22], the impact of DSTC in two-way 
amplify-and-forward relay channels was characterized. 

One of the principal challenges in distributed cooperative networks is to leverage the benefits 
from both NC and DSTC. A promising solution that addresses this challenge is space-time network 
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coding (STNC), which was proposed in [23]. Fundamentally, STNC combines the information from 
different sources at a relay, which involves the concept of NC. Moreover, STNC transmits the 
combined signals in several time slots using a set of relays, which involves the concept of DSTC. 
Based on these, STNC achieves spatial diversity with low transmission delay under the impact of 
imperfect frequency and timing synchronization. We note that in [23], the performance of STNC 
was evaluated for Rayleigh fading channels, where no closed-form expression was presented. 

In this paper, we consider a distributed cooperative network using STNC over independent but 
not necessarily identically distributed (i.n.i.d.) Nakagami-m fading channels, which generalizes the 
result in [23]. In this network, multiple sources communicate with a single destination with the 
assistance of multiple relays. Here, we focus on decode-and-forward (DF) protocol at the relays, 
which arises from the fact that this protocol has been successfully deployed in practical wireless 
standards, e.g., 3 GPP Long Term Evolution (LTE) and IEEE 802. 16j WiMAX [24]. Different from 
[23], we examine two fundamental questions as follows: "1) What is the impact of STNC on 
the symbol error rate (SER) in general Nakagami-m fading channels?" and "2 ) Can we provide 
closed-form expressions for the SER in Nakagami-m fading to alleviate the burden of Monte 
Carlo simulations ?" The rationale behind these questions is that Nakagami-m fading covers a 
wide range of typical fading scenarios in realistic wireless applications via the m parameter. 
Notably, Nakagami-m fading encompasses Rayleigh fading (m = 1) as a special case [25]. To 
tackle these questions, we first derive new closed-form expressions for the exact SER, which 
are valid for multiple phase shift-keying modulation (M-PSK) and M-ary quadrature amplitude 
modulation (M-QAM). To further provide valuable insights at high signal-to-noise ratios (SNRs), 
we derive new compact expressions for the asymptotic SER, from which the diversity gain is 
obtained. Specifically, it is demonstrated that the diversity order is determined by the Nakagami-m 
fading parameters of all the channels, but independent of the number of sources. Various numerical 
results are utilized to examine the impact of the number of relays, relay location, Nakagami-m 
fading parameters, power allocation, and nonorthogonal codes on the SER. Importantly, it is shown 
that nonorthogonal codes provide higher throughput than orthogonal codes, while guaranteeing full 
diversity over Nakagami-m fading channels. Our analytical expressions are substantiated via Monte 
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Carlo simulations. 



II. System Model 



Fig. 1 depicts a distributed cooperative network where L sources, Ui, U 2 , ■ ■ ■ ,Ul, transmit their 
own information to a common destination D with the aid of Q relays, R 1 , R 2 , . . . ,Rq. In this 
network, each node is equipped with a single antenna. We consider a practical and versatile 
operating scenario where the source-relay, the relay-destination, and the source-destination channels 
experience i.n.i.d. Nakagami-m fading. As such, we denote the Nakagami-m fading parameter 
between Ui and R q as mi q , the Nakagami-m fading parameter between R q and D as m qd , and the 
Nakagami-m fading parameter between Ui and D as m M . We further denote the channel coefficient 
between Ui and R q as hi q , the channel coefficient between R q and D as h q d, and the channel 
coefficient between Ui and D as hid, where 1 < I < L and 1 < q < Q. Throughout this paper, 
we define the variances of these channel coefficients as ~ CW (0, d^°), where G {l,q}, 
(p £ {q, d}, and ^ (p. Here, we incorporate the path loss in the signal propagation such that d^ 
denotes the distance between <fi and (p and a denotes the path loss exponent. 

In this network, the STNC transmission between the sources and the destination is divided into 
two consecutive phases: 1) source transmission phase and 2) relay transmission phase. In the source 
transmission phase, the sources transmit their symbols in the designated time slots. In this phase, 
the relays receive a set of overheard symbols from the sources. In the relay transmission phase, 
each relay encodes the set of overheard symbols to a single signal and then transmits it to the 
destination in its designated time slot. As illustrated in Fig. 2, (L + Q) time slots are required to 
complete the STNC transmission to eliminate the detrimental effects of imperfect synchronization 
on any point-to-point transmission in this network at any time slot. 

We proceed to detail the transmission in the two phases, as follows: 

In the source transmission phase, the signals received at the destination from Ui in the time slot 
I is given by 



where Pi denotes the transmit power at Ui, x t denotes the symbol transmitted by Ui, si (t) denotes the 




(1) 



spreading code of xi, and wi d is the additive white Gaussian noise (AWGN) with zero mean and the 
variance of N . The cross correlation between s p (t) and s q (t) are expressed as p pq = (s p (t) , s q (t)), 
where (/ (t) ,g(t)) = ^ f (t) g* (t) dt is the inner product between / (t) and g (t) during the 
symbol interval T. Moreover, we assume that p a = \\si (t)\\ 2 = 1. The signals received at R q from 
Ui is given by 

Vlq (t) = klqy/PlXlSl (t) + W iq (t) , (2) 

where w lq (t) is AWGN with zero mean and the variance of iV . 
In the relay transmission phase, the signal received at the destination from R q is given by 

L 

Vqd (t) = k qd ^2 Pql \fPq~i%l S l (t) +W qd (t) , (3) 



1=1 



where P q \ denotes the transmit power at R q and w q d is AWGN with zero mean and the variance 
of N . In (3), the scalar /3 q i denotes the state whether R q decodes xi correctly. Specifically, (5 q i is 
equal to 1 if R q decodes xi correctly, but otherwise. 

For the detection of the received signals at the destination, we assume that the full knowledge 
of the channel state information are available at the receivers with the aid of a preamble in the 
transmitted signal. We also assume that the destination has the detection states at the relays, which 
can be obtained via an indicator in the relaying signal. At the destination, the spreading codes s t is 
employed such that the information symbols x t is separated from y td and y qd , where I G {1, • • • , L}. 
For any desired symbol xi, the destination combines the information of x\ from Ui and the Q relays 
using maximum ratio combining (MRC). Therefore, the instantaneous signal-to-noise ratio (SNR) 
of xi is expressed as [23] 

^ = Pi \hldf + V- PqlPlqjhqdf ^ (4) 



q=l 



N ei 



where e\ is the Zth diagonal element of matrix R 1 associated with symbol x t and R is given by 



R = 



1 021 
Pl2 1 

PlQ P2Q 



PQ1 
PQ2 



(5) 



To facilitate the performance analysis in the following section, we re-express (4) as a unitary 
expression given by 

Q 

ll = C \h \ 2 + ^PqlCq \h q \ 2 , (6) 
9=1 

where c = Pid^ a /N denotes the equivalent SNR at D received from XJ h h denotes the unitary 
Nakagami-m fading coefficient between Ui and D with variance one, c q = P\d~^ /N ei denotes 
the gth equivalent SNR received at D from R q , and h q denotes the unitary Nakagami-m fading 
coefficient between R q and D with variance one. 

III. SER Analysis over Nakagami-tti Fading Channel 

In this section, we first derive new closed-form expressions for the exact SER with M-PSK and 
M-QAM. We then derive new compact expressions for the asymptotic SER, which will allow us 
to examine the network behavior in the high SNR regime. 

A. Exact SER 

In DF protocol, j3 q i denotes the decoding state at R q associated with x\. Based on the values of 
all P q i's, we define a decimal number as S t = [fin (5 2 i • • • /3qi] 2 to represent one of 2 Q network 
decoding states at Q relays associated with x t . Since all the channels in this network are mutually 
independent, the events that whether R q correctly decodes the received signal are independent. It 
follows that l3 q i's are independent Bernoulli random variables, the distribution of which is written 

1 - SER qh if (3 ql = 1 
G(M = { (7) 
j SER 9 /, otherwise, 

where SER g / denotes the SER of detecting xi at R q . Therefore, the joint probability of a particular 

combination of xi in Si is written as 

Q 

Pr(5,) = JlGCS„). (8) 

q=l 

Applying Bayesian rule, the SER of detecting xi at D is derived as 

2«-l 2«-l Q 

SER, = SER 7i|Si Pr (S{) = ^ SER 7!|S; JJ G ((3 ql ), (9) 

S t =0 S t =0 q=l 



where SER 7i|g; denotes the SER of detecting xi at D conditioned on Si. To facilitate the calculation 
of (9), we present the exact closed-form results for SER 7i|S; and G (/3 q i), as follows. 

1) Exact Results for SER 7l|s .- We commence the derivation of SER 7l|Si by presenting the PDF 
of 7j| Sj , / 7;|5( (v). If there are N "1" elements in one set Si, let a 1: a 2 , • • • ,a N G {ci, c 2 , • • • , cq} 
denote the equivalent SNRs of the N relays which decode xi successfully, where N < Q. We then 
define a = c to make the source equivalent to the zeroth relay. As such, the SNR at D to one 
desired symbol x t can be rewritten as 

N N 
2 ^ — 2 ^ — 2 

7i|Si = o \h \ + 2_^a n \h n \ = ) a n \h n \ . (10) 

n=l n=0 ^ 

J n 

According to [25], the PDF of F n in Nakagami-m fading is given by 

h - {y) = i^rw exp {~J ' 

where m n is the Nakagami-m fading parameter between the nth successful relay and the destination. 
For example, if the first and the second relays out of three relays successfully decode the information 
from Ui, we have mi = m w and m 2 —2d- Using Fourier transform together with [26, eq. (3.351.3)], 
the characteristic function (CF) of Y n is calculated as 

Cy n («) = r fy n (y) e^dy = ( 1 - J -^) . (12) 
Given that 7^ is the sum of Y n 's, the CF of 7^ is obtained as 

n=0 n=0 ^ m ™ / 

Applying inverse Fourier transform, the PDF of 7^ is derived as 



h lSl (*0 = 4 r ^ («) e-^d« 
z/l j —00 



00 



g( u ) 

We next seek the solution for g (u). Due to the randomicity of the wireless channels, we note that 
g(u) has N + 1 different poles z = —jm /a , z 1 = —jmi/ai, ■ ■ ■ , z N = —jm N /a N in complex 
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field. As such, based on residue theorem [27], the residue of k th pole of g(u) in complex field can 
be expressed as 

1 , d" 1 "- 1 



Res [g (z k ) , z k \ 



lim 



(m k — 1)! u ^ z k du mk 



- [(u - z k ) mk g («)] 



lim 



d m k -i 



{-j) m ka k m K (m k - 1)! u^z k du" 1 "- 1 



N 



,-juv 



n i- 



JUOn 

m n 



(15) 



As per the general Leibniz's rule, we derive Res [g (z k ) , z k \ in (15) as 

Res [g (z k ) , z k \ 



m k 



rn k 



mt-1 



i io l fc-2 *JV-2 



'm k - 1\ { i\ fi \ (i k -2\(ik-i 



io=Oii=0 «fc-i=0i fe+ i=0 ijv_i=0 



X 



X 



1 J Vw vi 

^«JV-2— ijv-i 



d u i N -2-iN-i 

(jik-2—ik-l 
d u i k -2-ik-i 



m N -i 



m k -i 



-mjv_i 



-m fc _i 



^*fc-l-*fe+i 



jua k+1 



du^- 1 %k + 1 \ in k+1 J 



-m k+1 



d io ~ h ( juaA 
du l °~ l1 \ mi J 



d 

du l 



\ io L _ jua Q \ 



Upon close observation, we simplify (16) as 

m k — 1 



Res [g (z k ) ,z k }= jB Njk:i v ' ' ' exp 



,m k -l-i 



i=0 



a k 



(16) 



(17) 



where B Nki is defined as 



B 



N,k,i 



m k mk (-l) 



l to l fc-2 »fc-l 



lN-2 



k V « y io=oil=0 j fc _ 1=0 j fc+1 =0 ijv-i=0 

\m k+ i) 



(18) 

i J VoJ \hj \ik-i) \ik+i) W-i 



x 



_ \ ijv-i / _ \ *JV-2— *JV~1 



m N 



J 



m N -i 



Qfc-i 



"ifc+l-*fc-i+«fe+l 



V m N a k J 



' t N-2~ l N 
-mjv-*JV-l 



xli r| ^fcQ/v-i \ " ' " " / 1 m k a k+1 \ 



m N „ia k 



-tk+l V" lv ±/ t k -2~ t k- 

-nijv-i — ijv-2+ijv-i 



X 1 



m fc _ia fe 



-™lfc-l-*fc-2+«fc-l 



... 1 



m x a k ) 



fnk+i^k J 

-mi— io+U / \ —mo— i+«o 

1-^ , (19) 



m a fe / 



and (m n ) i = Y (m n + i) /Y (m n ) is the Pochmann symbol. Based on the residues of the poles, we 
confirm that g(u) is available for the residue theorem. Specifically, using the residue of k th pole 
of g(u) in complex field, we obtain 

poo 

/ g (u)du = -2nj ^ Res \d ( z k) , z k ] . (20) 

The proof of (20) is shown in Appendix A. Based on (17) and (20), the PDF of 7i |s ; in (14) is 
derived as 

N m k -l , s 

A.., M=EE B N ^-" exp . (21) 

With the aid of / 7(|S; (u) in (21), we are capable to derive SER 7i|S; for M-PSK and M-QAM. 
First, SER 7i|S; for M-PSK is derived as 

I r(M-l)n/M roc / ^ 



SER 7!|giiM PSK = -y ^ a exp I - — ^ J / 7;|5( ( 7 )d 7 d0 

iV m fc -l f (M-l)n/M poo // h m , \ \ 



7T 

fc=0 i=0 

(22) 

where a and 6 are modulation specific constants. For M-PSK, a = 1 and 6 = sin 2 (7r/M). With 
the aid of [26, eq. (3.351.3)] and [28, eq. (5 A. 17)], (22) is derived as 

SER 7i| Si :MPSK 

a N mk ~ 1 r(M-l)it/M / ^ m \-m k +i 

= -Y] Y] B NAi Y{m k -i) / (— + — d6 

^ , n - n Jo \sm 2 9 a k J 



k=0 i=0 
N m fc -l 



iv m fc -± / x m fc -« j,^ -, -, / 7 / m k -i-i , x 

[^-WS£((^-o 5 ffl 

4(l + g) +sin(tan-^) £ (cob (tan" 1 a;))^ 1 



(23) 



where w = \ / °* & cot £ and T„ t = 



(?) 



We then derive SER 7!|S; for M-QAM as 

4 r' 2 r°° ( 6 7 



4 r /2 r°° / &7 \ 

SER 7;|SiiMQAM =- J y a exp I -— — Q 1 / 7i|5; (7)^ 

-H'T^i-^)^^ <24> 
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where a = 4 (l - 1/a/m) and 6 = 3/ (2 (M - 1)). Calculating the integrals in (24), SER 7;|S;i mqam 
is derived as 



SER 



7i| S( ,MQAM 
N rn k -l 



=4a J2 E S *.* r ( 



m k - i) 



k=0 i=0 



mi— % 



1 1 



, — + tan 1 u 1 ) > 

2 7T V a k b + m k \ \2 ' n 



IT 



2p 
P 



x 



4 ( 1 + ^jj +sin(tan-V) £ E T ^( 1 + ^j ( cos ( tan ~^)) 



1, , \\2(p-g)+l 



m, —1 



Aa2 H H B kji T(m k -i) I 

k=0 i=0 



0*\ 



m k 



1 1 



a k b 



m^—i—l 



4 71 V Ofcfr + \ ^ 2 



(I - tan ^2) 



p=0 



2p 
P 



x 4 1 + 



m k 



m—i—l p 



— sin (tan 1 w 2 ) ^^T P:q (l-\ — — ) (cos (tan 1 u}2)) 

p=l q=l ^ 



-p 



1 _ nn2( P - 9 ) + 1 



(25) 



where Ul = J tt cot ? and w 2 = J ?* b cot f . 

2) Exact Results for G(f3 g i): We now analyze G(/^) for M-PSK and M-QAM, respectively. 
According to (7), it is equivalent to analyze SER^. Using (2), the received SNR at R q is written 



as 



Pi IVI 2 ^gi" it, |2 1, |2 



(26) 



where = Pid ql a /N denotes the equivalent SNR at i? g received from Ui, and ft, denotes the 
unitary Nakagami-m fading coefficients between JJ\ and R q with variance one. 
We first derive SER,, for M-PSK as 



SER, 



<jZ,MPSK 



a(M - 1) a c q ib 



M 



sin tan" 



7r y c g /6 + miq 



+ tan' 



p=0 



2p ) f 4 fi + a» 

V ) \ \ m q 



mi q -l p . L\ -P 

1 y y T p t ( 1 + ^ ) (cos (tan- 1 w) ) 2(p -' )+1 



(27) 



where w 



cot 



c q ib+m lq M 
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We then derive SER qi for M-QAM as 

4a 



SER, 



<ji,MQAM 



~-2a - 



Cqlb 



7T \l c q ib + m k 

m~i—l p 



rrik— i— 1 



(^ + tan V) Yl 



p=0 



2p ) U(i + C st 

Pj\\ rn k 



+ sin (tan-V) ^ 1 + -^ (cos (tan" 1 ^)) 

p=l q=l V mfe / 



-a 2 + 



4a 2 



a k b 



7r \ a k b + m k 

m—i—l p 



rrik—i—l 



(| - taI1 ^2) ^ 



p=0 



2p 
V 



m-i-i p , c b\~ P 

— sin (tan _1 n7 2 ) ^ ^ T M ( 1 + j (cos (tan _1 tJ7 2 )) 

p=l 9=1 'V 



2(p-?)+l 
1 + ^ 

m k 

2(p-?)+l 



(28) 



where 



a k b+m u 



cot ^ and ro 2 



-g^cotf. 

a k b+m lq 4 



Substituting (27) and (28) into (7), we obtain G for M-PSK and M-QAM, respectively. 
Therefore, we insert (7) and (23) into (9), which yields the exact closed-form SER for M-PSK, 
and substitute (7) and (25) into (9), which gives the exact SER for M-QAM. Observing (23), 
(25), (27), and (28), we see that the exact SER expressions for M-PSK and M-QAM are given in 
closed-form and are valid to arbitrary numbers of sources and relays. 



B. Asymptotic SER 

We now provide useful insights into the network behavior in the high SNR regime. In doing so, 
new compact closed-form expressions are presented for the asymptotic SER. 
We first focus on M-PSK. Based on (10), SER 7;|S; for M-PSK can be alternatively written as 



SER 7(| mpsk 

r(M-l)w/M poo 



a 
k Jo 
a f 
n Jo 



o exp {-^e) ^ ^ d ^ dd 



(M—l)ir/M I poo poo 
JO 



N 



N 



, exp \- bJ Bf ) (n^^)jn^)^ (29) 



Substituting (11) into (29) and using [26, eq. (3.351.3)], we obtain SER 7!|S iMPSK as 



SER 



-7; | S , ,MPSK 



f (M-l)n/MN 

71 J* to 

r (M-l)n/M N 
= ~ f 

71 J° to 



J a^Y (mi) 



exp 



+ — ) Vi ) dy, t ) d9 



sin 9 cii 



a™ % Vsin 2 ^ ' at 



1 + Ue. 



(30) 
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We next use [26, eq. (2.513.1)] to develop an asymptotic expression for SER 7;|5 )M psk 

SER 7 - Si , MPS K n fyr (-^) ) dd 



as 



vi=0 

AT / \ m; 



i=0 



= -^5— ^n(-) ' < 31 > 



where 



1 ^Ef = 0^ (M-l)7T (-igj^ 

^ " 2 2 E£o - V Ef=o M + ytfU 

I . — )\r 



x E (-D* H°" T U „ LJLJ - <32) 

Similarly, an asymptotic expression for SER^mpsk is obtained as 

opRoo ^ a (m lq yn lq A M:miqfi 

^^,MPSK S - J m lq ■ 

Correspondingly, the asymptotic G (/3 qt ) for M-PSK is given as 

I 1, if Pal = 1 

G (^)mpsk = (34) 
I sef V,mpsk> otherwise. 

Based on (31), (33), and (34), the asymptotic SER for M-PSK is derived as 

2«-l Q 

SER-^mpsk — ^ SER^ ff _ >MPSK ]~[ G (^)mpsk • (35) 

5 i= 9=1 

Following the same procedure outlined for M-PSK, we derive the asymptotic SER for M-QAM 

as 

2«-l Q 

SER^qam = X] ^ ER 7^5 ; ,MQAM II & (^)mQAM ' ( 36 ) 
S ; =0 9 =1 

where the asymptotic SER 7;|5;i for M-QAM is derived as 

SER^mqam < z^^" Awf[ --^—A 4/3 , N f[(^j , (37) 

i=0 x ' 1=0 x 7 

with A 2: mi,N and Ai/3 )J7lj) jv being defined in (32), and the asymptotic G (f3 q i) for M-QAM is derived 



as 



,1, if Pgl = 1 

G (AzOmqam = ^ ( 38 ) 
SER oJ,mqam> otherwise, 
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with 

QffB co ^ 4a fm lq \rn lq A2^fi W (mi q yn lq A A/3>miqfi 
b ^iMQAM<—{—) -nzr —{—) c m q ■ W 

Based on (35) and (36), we next examine the diversity order of the network, which represents 
the slope of the SER against average SNR in a log-log scale. According to (35) and (36), the 
asymptotic SER of Ui can rewritten as 

2Q-1 

sERr = e ^ 

Si=0 



N 1 Q 



q=i c qi 



J=0 

where Q Sl denotes the coefficient independent of a, L and c^.in which the diversity order can be 
confirmed as 

div = m ld +) ^ min (m qd , m lq ). (41) 

It is evident from (41) that the full diversity order is achieved, which is determined by the Nakagami- 
m fading parameters of all the channels. Notably, the diversity order is independent of the number 
of sources. In particular, this full diversity order is preserved even non-orthogonal STNC codes are 
employed. 

IV. Simulation and Numerical Results 

In this section, simulation and numerical results are presented to examine the impact of network 
parameters with STNC (e.g., the number of relays, the relay location, Nakagami-m fading param- 
eters, and power allocation) on the SER of Ui. In the figures, we consider a practical scenario 
where the relays are placed at different distances from D and Ui with cj ^ c.i and Cji ^ c a for 
j i. We set the distance between C/j and D as di d = 1. The cross correlations between different 
spread codes, defined in (1), are set to be zero. We also assume equal transmit power at each node. 
Further, our results concentrate on the practical example of a highly shadowed area with the path 
loss exponent as a = 3.5 [29]. In the figures, the exact SER for M-PSK is evaluated by substituting 
(7), (23), and (27) into (9), and the exact SER for M-QAM is evaluated by substituting (7), (25), 
and (28) into (9). The asymptotic SER for M-PSK and M-QAM is calculated from (35) and (36), 
respectively. 
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A. Impact of Number of Relays and Equal Nakagami-m Fading Parameters 

In this subsection, we focus on equal Nakagami-m fading parameters with m 8 = m. The average 
received SNRs are set as C( i+ i)i = c i+ i = Ci + O.I7A and c = c z = O.67A. Fig. 3 plots the exact 
and asymptotic SER with 4QAM. Fig. 4 plots the exact and asymptotic SER with 8PSK. From 
Figs. 3 and 4, we see that the asymptotic SER curves accurately predict the exact ones in the high 
SNR regime. By observing these asymptotic curves, it is evident that the diversity order increases 
with Q, which indicates that increasing the number of relays brings an improved performance. It 
is also seen that the diversity order increases with m, which indicates that the improvement in 
fading channels leads to a reduction in the SER. Moreover, we see that the simulation points are 
in precise agreement with our exact analytical curves, which demonstrates the correctness of our 
analysis in Section III. Comparing the SER in Fig. 3 with that in Fig. 4, we further see a poorer 
network performance is achieved by higher order modulation schemes. 

B. Impact of Relay Location 

In this subsection, we consider d tq 7^ d qd , which leads to c q 7^ c qd , and consider equal Nakagami- 
m fading parameters with m 8 = m = 2. We further normalize d td to unity with d td = 1. Fig. 5 
plots the exact SER with BPSK for Q — 2. In this figure, Cases 1, 2, 3 represent the scenario 
where the relays are located close to the source, while Cases 4, 5, 6 represent the scenario where 
the relays are located close to the destination. 

We first consider Cases I, 2, and 3. We see that Case I offers a prominent SNR advantage relative 
to Case 2. This indicates that the reduction in the distance between the relay and the destination 
brings a substantial SER improvement. We also see that Case I and Case 3 achieve almost the 
same SER across the entire SNR range. This indicates that the SER improvement from the reduced 
distance between the source and the relay is negligible. These observations are due to the fact that 
the network performance is dominant by the relay-destination link when the relays are close to the 
source. As such, the quality improvement of the relay-destination link has a higher positive impact 
on the SER than that of the source-relay link. 

We next consider Cases 4, 5, and 6. It is seen that Case 4 provides a substantial SNR advantage 
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compared to Case 5. It is also seen that Case 4 achieves a slight SNR advantage compared to Case 
6. These observations are explained by the fact that the network performance is dominant by the 
source-relay link when the relays are close to the destination. 

C. Impact of Unequal Nakagami-m Fading Parameters 

We concentrate on unequal Nakagami-m fading parameters and set the average received SNRs 
as C( i+ i); = Cj + i = Ci + O.I7A and c = c ; = O.67A. Fig. 6 plots the exact SER with 4QAM for 
Q = 2. This figure clearly shows that the diversity order in (41) is accurate. For example, it is 
evident that the asymptotic SER curves of Cases 1, 2, and 3 are in parallel, which indicates that 
they achieve the same diversity order. As indicated in (41), Cases I, 2, and 3 achieve identical 
diversity order of 3. Moreover, we see that the diversity order of Case 4 increases to 4, and the 
diversity order of Case 5 increases to 6. This is predicted by (41), which shows that the diversity 
order is determined by the Nakagami-m fading parameters of all the channels. 

D. Impact of Power Allocation 

We now focus on arbitrary transmit power at each node. We consider equal Nakagami-m fading 
parameters with m 8 = m = 2, set the relay location as d tl = 0.8, d t2 = 1, d ld = 0.9, and d 2d = 0.7, 
and normalize d ld as d td — 1. We denote the transmit powers at C/j, R±, and R 2 as P , Pi, and P 2 , 
respectively. Under the total power constraint, we have P + Pi + P 2 = 3P. Fig. 7 plots the exact 
SER versus f = Pi /(Pi + P 2 ) with 4QAM for Q = 2 and P/N = 12 dB. We see that the optimal 
value of P\jP 2 depends on P . For example, the optimal power allocation is at Pi = 1.75P and 
P 2 = 0.75P when P = 0.5P. Moreover, the optimal power allocation is at Pi = P 2 = 0.5P when 
P = 2P. Using our SER expressions with different relay locations, iterative search method can 
be used to find the optimal power allocation that minimizes the SER. 

E. Impact of Nonorthogonal Codes 

We now turn our attention to the impact of nonorthogonal codes. We consider p pq = p ^ 
for all p, q and equal Nakagami-m fading parameters with m; = 2. Fig. 8 plots the exact SER 
with 4QAM for Q = 2 and N — 3. The case of p = represents orthogonal codes. We see a 
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reduction in the SER as p increases. We also see that the diversity order is not affected by cross 
correlation. As such, the nonorthogonal codes which permit broader applications can be used for 
higher throughput without sacrificing the error rate significantly. 



In this paper, we analyzed the SER of STNC in a distributed cooperative network where L sources 
communicate with a single destination with the assistance of Q relays. For M-PSK and M-QAM 
modulation, new exact closed-form expressions of SER over independent but not necessarily iden- 
tically distributed Nakagami-m fading channels were derived. Moreover, the asymptotic SER was 
derived to reveal the network performance in the high SNR regime. Specifically, the asymptotic SER 
reveals that the diversity order of STNC was determined by the Nakagami-m fading parameters of 
all the channels. Simulation results were used to validate our analytical expressions and to examine 
the impact of Nakagami-m fading parameters, relay location, power allocation, and nonorthogonal 
codes on the SER. 



According to Fig. 9, all poles are located in a closed curve. According to the residue theorem, 
the integral over a closed curve can be expressed as the linear combination of the residues of the 
poles in the curve. Mathematically, we have 



V. Conclusions 



Appendix A 



Proof of (20) 




(42) 



where Cr is the counterclockwise semicircular curve in Fig. 3. The absolute value of J Cr g (z) dz 
in (42) can be upper bounded as 




17 



r N 



dz. 



(43) 
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When R — > oo, which indicates that the integral range of g (u) is from oo to — oo, we find the 
property of (43) as 
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(44) 



Substituting (44) into (42), we obtain (20) and thus complete this proof. 
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Fig. 1. System model. 
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Fig. 2. A framework of space-time network coding. 
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Fig. 4. Exact and asymptotic SER with 8PSK for Q = 1, 2, 3, cy+i); = Ci+i = a + O.I7A, and co = Co; = O.67A. 



21 





Fig. 6. Exact SER with 4QAM for Q = 2, C( i+ i); = c»+i = d + O.I7A, co = coi 
Case 2: run — l,mi2 = l,mid = l,m\d = 2, m,2d = 2; Case 3: run = 2, mz2 
ma = 1, m;2 = 2, mid = 1, mid = 1, m-id = 2; Case J: mj = m = 2. 



= O.67A and 6 cases: Case 2: mi = m = 1; 
= 2, m;d = l,mid = l,m2d = 1; Case 4: 
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Fig. 9. The distribution of the poles. 



